The mechanism of phase transitions of infinite-range systems enlightened by an 
elementary Z 2 -symmetric classical spin model 
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In this paper we present the most elementary model that we know with a first order Z2-symmetry 
breaking phase transition. It is a classical spin model with potential energy assuming only two 
values, which despite its dramatic simplicity reproduces all the characteristic features of a first order 
symmetry breaking phase transition. Its aim is not to reproduce anyway some physical systems, but 
to enlighten the mechanics at the deep origin of a phase transition of any system. We consider the 
model as elementary in this sense. Indeed, it reveals in the most evident way how the Boltzmann 
factor competes with the entropic factor in order to generate the phase transition by varying the 
temperature. The paper ends revisiting the solutions of the Ising model and the spherical model 
(Berlin-Kac) in the mean-field version which show the same picture of the model introduced in 
this paper, but extended to continuous phase transitions. A limit of this analysis is that all the 
three models here considered satisfies the strong constraint of having the potential expressible as a 
function of the magnetization, thus they belong to the class of infinite-range systems. Anyway, we 
hope that this study may be helpful in finding out the general sufficiency conditions under which a 
potential can entail a phase transition in a general physical system also in the finite-range case. 

PACS numbers: 75.10.Hk, 02.40.-k, 05.70.Fh, 64.60.Cn 
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I. INTRODUCTION 

Phase transitions are sudden changes of the macro- 
scopic behavior of a physical system composed by many 
interacting parts occurring while an external parameter is 
smoothly varied, generally the temperature, but e.g. in a 
quantum phase transition is the external magnetic field. 
The successful description of phase transitions starting 
from the properties of the microscopic interactions be- 
tween the components of the system is one of the major 
achievements of equilibrium statistical mechanics. 

From a statistical-mechanical point of view, in the 
canonical ensemble, a phase transition occurs at spe- 
cial values of the temperature T called transition points, 
where thermodynamic quantities such as pressure, mag- 
netization, or heat capacity, are non-analytic functions 
of T. These points arc the boundaries between different 
phases of the system. Starting from the exact solution 
of the two-dimensional Ising model [l7[ by Onsager [27| . 
these singularities have been found in many other models, 
and later developments like the renormalization group 
theory [l5| have considerably deepened our knowledge of 
the properties of the transition points. 

But in spite to the success of equilibrium statistical me- 
chanics the issue of the deep origin of a phase transition 
remains open, and this motivates further studies of phase 
transitions. In this paper we present a Z2-symmetric clas- 
sical spin model showing a first order symmetry break- 
ing phase transition which, in our knowledge, is one of 
the most elementary models capable to show such a phe- 
nomenon in the canonical treatment. We think that in 
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it the general mechanism of entailing symmetry break- 
ing phase transitions acts, and that it becomes manifest 
due to the extreme simplicity of the model. The model 
is derived from the " hypercubic model" introduced in Q 
by replacing the real variables with classical spin ones, 
and it can be consider as a drastic simplification of the 
mean-field Ising model. 

In more physically realistic models we think that the 
mechanism entailing phase transitions may be the same, 
but that it remains hidden due to the complexity of the 
models, also in that cases in which the canonical ther- 
modynamic is analytically solvable. Indeed, even if the 
analytic solution is at disposal, we can find out whether 
a system undergoes to a phase transition or not only by 
analyzing the solution, but we do not know any general 
criterion capable to predict such a phenomenon founded 
only on the properties of the potential energy landscape. 

Anyway we have to remark a limitation of the model 
we present, because being the potential writable as a 
function of the order parameter, a.i. the magnetization, 
the model has to be considered belonging to the class of 
infinite-range systems, which are not very realistic from 
the physical viewpoint, mostly in classical physics. Any- 
way, the mechanism entailing the phase transition may 
contain some features which may be in future generalized 
to the finite-range case. This is under study. 

Returning to the potential, it is a double- well function 
of the magnetization m. Since the entropy s(m) is a con- 
cave function due to the geometry of configuration space, 
the phase transition emerges by the competition between 
the potential and the entropy in shaping the free energy 
f(m,T) — v — Ts, which can have one or two absolute 
minima with respect to the magnetization depending on 
the value of the temperature T. If T is less than the 
critical temperature, / becomes a double-well function 
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as well the potential and the symmetry is broken, while 
if T is greater than the critical temperature, / remains 
a convex function as well the entropic term —Ts and the 
symmetry remains unbroken. 

This picture of the free e nerg y is exactly the one imag- 
ined by Landau-Ginzburg jl5l fioj ] in order to describe 
phenomenologically the classical continuous phase tran- 
sitions, but it can apply also to the hrst order ones as for 
the model of this paper. 

Section [IT] is devoted to the introduction of the new 
model and to its study in great details. Further, Section 
IIII1 and IIVI are devoted to the Ising pjj and Berlin-Kac 
[3| models in the mean-field version revisited in order to 
show the occurrence of their well known classic contin- 
uous phase transitions by the same mechanism of the 
model introduced here. 



II. THE MODEL 

In a recent paper Q a simple model with first order 
Z2-symmetry breaking phase transition has been intro- 
duced. It has been called "hypercubic model" because 
the interaction among the N continuous degrees of free- 
dom is by a double- well potential with a gap proportional 
to N constructed by hypercubes of M . 

From the hypercubic model we derive a more essential 
model by replacing the ith coordinate qi with the classical 
spin <Ji 6 {— 1,+1}, the same of the Ising model. In 
order to reproduce the double-well potential we define 
the potential in the following way 



V(a) = 



—NJ if (Tj = 1 and a% 
otherwise. 



-1 Mi 



(1) 



In this way the potential has a gap proportional to N 
which separates the two wells. The constant J plays the 
role of the strength of the interaction, and it modulates 
the critical temperature of the phase transition. 

From a physical viewpoint such a potential may be re- 
garded as the tendency of the spins to be in the two states 
{1, ■ ■ • , 1}, { — 1, • • ■ , —1} by a completely delocalized in- 
teraction among the spins proportional to their number 
N . This is similar to what happens in some quantum 
systems, where the non-locality of the wave function can 
produce particular kinds of interaction depending only 
on the number of particles and not on their distances. 

The thermodynamic can be solved in one shot, by using 
the decomposition formula 




FIG. 1: From topo to bottom, and from left to right. Free 
energy, average potential, specific heat, and magnetization of 
the model introduced in section[TT]for J = 1. The smooth lines 
are for TV = 5, 10, 30, and show the nonuniform convergence 
toward a discontinuous limit corresponding to a first order 
phase transition with T c — 



In the limit of large N only one between the two exponen- 
tials of the right hand side of the last equation survives, 
so Zpf can be approximated as 



Zn ~ * if T <T P 



(4) 



where T c — is the critical temperature. T c separates 
two different analytic forms of Zn(T). 

In the thermodynamic limit the free energy, the aver- 
age potential, and the specific heat are respectively 



v = -T' 



d (I 



dT \T 



dv 
dT 



-J if T < T c 
-Tin 2 if T > T r 



if T > T c 
-J if T < T c 



if T > T c 
+oo if T = T c 
if T < T c 



(5) 



(6) 



(7) 



Z N =J2 e ~ Z ^ L =X)e-£w JV (K) (2) 
{*} v, 

where uj^(V) is the microcanonical density of states at 
fixed N, so 



We can also give the relative expressions for finite N 
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In ( 2^ - 2 + 2e 



N 



-) (8) 



Z N = e-Tuj(0) + e—Lu(-NJ) 
= 2 N -2 + 2e^ =e Nln2 - 



N J 

2e~. 



(3) 
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-2Je t 



2 N -2 + 2e^ 



(9) 
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C V ,N 



dvN 
dT 



2J 2 Ne^ (2 N -2) 



2 N -2 



2e t 



rp2 



(10) 



Now consider the magnetization. As T < T c the aver- 
age potential is v = — J, so that the only configurations 
accessible to the representative point are {1, • • • ,1}, and 
{— 1, • • • ,—1}, at which the values of the magnetization 
per degree of freedom m = 1, and m = — 1 correspond, 
respectively. Since in the thermodynamic limit the prob- 
ability of overturning simultaneously all the spins is van- 
ishing, the magnetization has to be fixed in one of the 
two possible values. 

As T > T c the average potential is v = 0, so that 
the representative point can freely go around the whole 
lattice expect the sites {1, • • • ,1}, and {— 1, ••• ,—1}. 
The magnetization is vanishing because in the thermo- 
dynamic limit the most probable configuration is the one 
with half spins 1, and half spins —1. Thus 



±1 




if T <T C 
if T > T c 



(11) 



Summarizing, we are in front of the complete picture 
of a first order ^-symmetry breaking phase transition, 
despite the dramatic simplicity of the model. 

It is worth noticing that the potential ([1} has to be 
considered as the effect of a completely delocalized inter- 
action among the spins, thus it is infinite-range. 



has to hold between the volume of the hypercubes A + , 
or A~ , and the hypercube B, so that bja — 2 has to be 
assumed. Finally, in order to complete the mapping we 
make the identification v c i— > J. 

The critical temperature results T c = v c / \n(b/a) > 
1/ In 2, as it has to be, and the thermodynamic functions 
are the same. 

Since the model of this section depends on one param- 
eter only, it is in some extents more essential in what 
concern the mechanism giving rise to the phase transi- 
tion with respect to the hypercubic model. Further, due 
to the discreteness of configuration space, it is possible 
to evaluate the effect of an external magnetic field, as we 
will see in Section Hi D[ 



B. The free energy as a function of m and T 

In Section [TT] we have calculated the partition function 
by decomposing it in a sum on the values of the potential 
which are 0, —NJ, but we can also choose to decompose 
it in a sum on the values of the magnetization, which will 
begin an integral in the limit N — > oo. In this way we 
can deduce the shape of the free energy with its minima 
as a function of the magnetization, and moreover we can 
evaluate the effect of an external magnetic field applied 
to the system. 

Let k be a positive integer which counts the number 
of spins <7i = —1 in a configuration of the system, so 
< k < N. k labels the subsets of configuration space 
{ — 1, 1} N at constant magnetization defined as 



A. Mapping in the class of the hypercubic models 

We recall that the configuration space of the hyper- 
cubic model presented in [2[ is made by a hypercube B 
centered in the origin of coordinates where the poten- 
tial takes the value v c > 0, except in two smaller hy- 
percubes A + , and A~ included in B where the potential 
takes the value 0. The hypercubes A + , and A" are such 
to be the images one of the other under Z2-symmetry. 
The model shows a first order symmetry breaking with 
T c = v c / ln(6/a), where b is the side of the hypercube B, 
and a is the side of the hypercubes A + , and A~ . So, the 
class of the hypercubic models depends on two free pa- 
rameters: v c € [0, +oo), and the quotient b/a £ [2, +oo). 
If b is strictly greater than 2a, some freedom to dispose 
A + , and A~ into B advances, but it does not affect the 
thermodynamic functions, but only the value of the mag- 
netization. 

The class of the models introduced in this section de- 
pends on the parameter J 6 [0, +00). The mapping from 
this class and the hypercubic models one can be made by 
identifying the lattice sites {1, • • • ,1}, and {— 1, • • • , — 1} 
with the centers of mass of the hypercubes A + , and A~ . 
Since the statistical weight of a single lattice site is 1 and 
the weight of the whole lattice is 2^ , the same proportion 



{-1,1} 




(12) 



The configuration space is reconstructed by the disjointed 



union of all the Sfc's: {—1, 1}^ = Ufclo The relation 
which links k to the corresponding magnetization m/j is 
rnfc = 2k/N — 1, so that —1 < rrik < 1- 

The potential can by written as a function of k 



-NJ if k = 0,N 
if < fc < N 



(13) 



and so the partition function can be decomposed as a 
summation over all the values of k 



N 



z N = Y J < 



-pv k 



vol(Ek) 



(14) 



k=0 



where vol(Ek) ("vol" stand for "volume", even though 
this word generally refers to continuous sets, but with a 
small abuse of language we use it also for discrete sets. A 
more appropriate word may be " cardinality" ) is the num- 
ber of lattice sites belonging to vol(Yik) is nothing 
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but the density of states at fixed magnetization, which is 
liked to the entropy density at fixed magnetization Sk by 
the relation uoZ(£/-) = e * . 
It turns out that 



voZ(E fc ) = 



N 



k\(N-k)\ 



where (,) are the binomial coefficients. So 



(15) 




FIG. 2: Free energy as a function of m as T < T c (left), and 
as T > T c (right). 



E 



N-l 



,N/3J 



E 



Z N = > i 

k=0,N k=l 
N-l 



Nl 



kl(N-k)\ 
2e ^+^e ln ^™. (16) 



fe=X 



In order to search for the thermodynamic limit, we can 
use the Stirling formula to approximate the factorials in 
the binomial coefficients for large N 



N\ 




( (2tt) 3 N 



kl(N-k)\ \k{N-k) 



Recalling that j| = mfc 2 +1 , we can express Zjq in terms 
of the magnetization rrik , and then make the substitution 
TOfc i y to where to is a continuous real variable belong- 
ing to the interval [—1,1]. This substitution is possible 
because as N — > oo the values of to* becomes dense and 
equally spaced in the interval [—1, 1]. Thus 



E 



,NPJ 



m=-l,+l 



dme~' 9y(m) woZ(S„ l ) = 



Writer 

am e \ 



(18) 



The potential, and the entropy densities as functions of 
to remain defined respectively as follows 



— J if to = ±1 
if - 1 < to < 



(19) 



Remark. The fact that the potential (JlJ is a function 
of to lets to separate the contribute of the Boltzmann fac- 
tor and the density of states at fixed to in the calculation 
of the partition function. Since the density of states in- 
creases as e N due to the dimension of configuration space 
which increases as N, setting the jump proportional to N 
is a necessary and sufficient condition in order to entail 
the phase transition. Indeed, if the gap increased faster 
than N the potential would not be intensive, while if the 
gap increased more slowly than N the phase transition 
would disappear. 



C. Fisher zeros 

The fact that the partition function can be exactly 
solved allows us to discuss the exact localization of the 
Fisher zeros [l2|], the zeros of the partition function in 
the complex temperature plane. Because the analyticity 
of Zn, the zeros have not to lie on the real axis, but if 
a singularity is present in Z in the limit ./V — > oo, they 
have to converge to the real axis in correspondence of the 
critical temperature. The Fisher zeros are the analog of 
the zeros of the grand canonical partition function in the 
complex fugacity plane introduced by Yang and Lee in 
25]. 

We introduce the inverse temperature j3 = l/T and 
solve the equation 



Z N (P) = 



(21) 



whose solutions are 



, = ,»( » (i^)'"V (20, 

\ 1 — TO \ 1 + TO / / 

The free energy / = v — Ts is plotted in Figure [5] 

The central minimum of / is —Tin 2 which competes 
with the two other extreme minima — J in order to de- 
termine the absolute minima, and so giving rise to the 
critical temperature of the system T c — ■ The mag- 
netization is already given in (|lip . 



= i(2k + 1)tt 1 , N _i _ 
ho NJ NJ K V ' 

As we expected, the limit TV — > oo of each solution is the 
inverse critical temperature 



Jim 03 o ) = /?c=^. (23) 
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D. Effect of an external magnetic field 

In this section we determine the effect of an external 
magnetic field applied to the system. The new potential 
has to take into account the magnetic interaction by the 
term 




JY 



°%H = -NmH 



FIG. 3: Effect of an external field on the free energy (|26|) as 
(24) H > (left), and as H < (right). The temperature is below 



so the new potential is 

V H =Nv- NmH (25) 
where v is given in (|19[) . The new free energy is therefore 



fu = v — Ts — mH 



(26) 



where s is given in (1201) . fjj is plotted in Figure [3l 

In order to find out the magnetization we have to min- 
imize fu with respect to m. We start by solving the 
following equation 



df H _ 
dm 

whose solution is 



T 



In 



H = 



(27) 



m cent = tanh [ — 



(28) 



Then we confront /ii (m ce „ t ) with /fj (1) if H > 0, or with 
/h (— 1) if H < 0. This can be easily made by introducing 
the quantity 



m 




FIG. 4: Magnetization as a function of T for some values of 
the external field H = 0,0.5, 1,2 (red, green, blue, magenta). 
J = 1 (continuous lines), and J = (dashed lines). In the 
inset the function (|29l) . whose root is the pseudo-critical tem- 
perature T C (J,H) which separates the two analytic forms of 
the magnetization. 



III. THE MEAN-FIELD ISING MODEL 



A/ff = /ff(mcent) — /h(±1) 

= -T c In (l + e 2 ± (2H + J) = 0. (29) 

respectively if H J* 0. The root gives a sort of pseudo- 
critical temperature T c ( J, 7?) which separates two differ- 
ent analytic forms of the magnetization 

_/ ±1 if T<T C (J,H) 

m -\±tanh(§) if T>T C (J,H). W 

In the limit H i— > ±oo, T C (J, H) > oo, so the magneti- 
zation remains frozen in m = 1 or m = -1 for all values 
of T. 

We observe that the presence of the external field 
breaks the symmetry of the system, but does not elim- 
inate the phase transition meant as the presence of a 
nonanalyticity in the thermodynamic function. This is 
not surprising, and it is due to the fact that the phase 
transition is of the first order, and not continuous. 



The solution of the classical mean-field Ising model is 
by means of the mean-field theory and it is a paradig- 
matic example in literature, e.g. Nevertheless, 
we revisited it enlightening how it lets to separate the 
contribution of the potential from the entropy in the free 
energy at fixed magnetization, showing as they compete 
in order to entail the classical second order phase transi- 
tion occurring in this model. 

The potential is as follows 



v = -j^T, a ^ ( 31 ) 

where cr; e { — 1,1}, J > 0, and the factor 1/N is intro- 
duced in order to guarantee the intensivity of the poten- 
tial per degree of freedom. 

Following the same method already used to solve the 
model in Section [TTl we introduce the positive integer 
k which labels the subsets of the configuration space 
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FIG. 5: Free energy as a function of the magnetization of the 
mean-field Ising model (J3SJ as T < T c (left), and as T > T c 
(right). 



{ — 1,1}^ at constant magnetization m k whose elements 
have k positive spins. V can be easily written as a func- 
tion of k because the interaction among the spins is "all 
with all" 

V k = --(N-2kf = -JN[l--j . (32) 

Since k is linked to the magnetization by the relation 
TOfc = 2k /N — 1, the potential can be expressed as a 
function of m k 



IV. THE MEAN-FIELD SPHERICAL MODEL 

In this Section we add another example of a mean- 
field model solvable in the same way applied in previous 
Sections. It is the model introduced by Berlin and Kac 
[3, HH , and also known as the spherical model. It approx- 
imates the feature of the Ising model by substituting the 
discrete classical spin variables Oj with the continuous 
real qi, and constrained them on a iV-spherc of radius 
y/N, where N is the number od degrees of freedom. Such 
an TV-sphere contains the lattice sites of the Ising model 
if a t £ {-1,1}- 

The potential is the same of the Ising model 



v = -JjJ2 Ji 39iQj ( 37 ) 

<>3 =1 

where the sum is extended up all the couples of variables 
because we consider only the ferromagnetic mean-field 
version, and so we also set J^- — J > 0. The factor l/N 
has been introduced in order to guarantee the intensivity 
of the potential per degree of freedom. 

Thanks to this choice, V can be written as a function 
of the magnetization m — l/^X^Li Qi 



V k = -JNmi 



(33) 



Then, as N 
where to G 
finally 



-> co we can made the substitution mj. — > m 
—ljl] is a continuous real variable, getting 



V 

N 



-Jm . 



(34) 



The entropy is the same of the model in Sectionllll already 
found out in (|20[) . So the free energy results to be 



/ = -Jm 2 - Tin 



1 



1 — to 
1 + to 



(35) 



/ is plotted in Figure [5] 

In order to find out the spontaneous magnetization we 
have to differentiate / with respect to m and sets to zero, 
the resulting equation gives m(T) in the implicit form 



= tanh 



2Jm 



(36) 



The last equation has two symmetric solutions as T < 2 J, 
and as T > 2 J, so the critical temperature is T c = 2 J, 
as well known. 

We do not add nothing about the effect of an external 
filed because it is similar to what happens in the case 
of model in SEC. II, and because it is well known in 
literature, e.g. [16j |. 




= -Jm 1 



(38) 



and so the free energy as a function of m and T can 
be derived. In order to do this, we want to calculate 
which is the intersection between the 
= Nm and the iV-sphere of radius 
N. Thus, S m is an (N — l)-sphere and its volume is 



the measure of S m 
hyperplane Y^iLi * 



vol(E m ) = 



2n~ 



r(m) 



N-2 



(39) 



where the radius r(m) is a function of the magnetization. 
r, by the Pythagorean theorem, is 



(40) 



r(m) = VN(1 - tti ) - 



where m belongs to the interval [—1,1]. The sketch in 
Figure [S] can help the comprehension. 
The entropy results to be 



sn 



1 

N 



hi vol (S m ) 



(41) 



and by applying the Stirling approximation of N\ we ob- 
tain in the limit N — > oo 



s(m) = - ln(l 



to 2 ) + — In - 



(42) 
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FIG. 6: Mean-field spherical model. The circle represents the 
configuration space schematically drown at N = 2, and the 
secant (red) represents the E m whose radius r is linked to the 
magnetization m via the Pythagorean theorem applied to the 
triangle sketched. 




m 



FIG. 7: Free energy as a function of the magnetization of 
the mean- field spherical model (|43[) as T < T c (left), and as 
T>T C (right). 



Finally the free energy f = v — Ts results 

/ = - Jto 2 - j ln(l -m 2 )-j In ^PhTe. (43) 
By minimizing with respect to to we find out m(T) 

m ={ ± ( l -hY if T < T c (44) 
[0 if T > T c 

where T c = 2J is the critical temperature. rn(T) shows 
the well known continuous phase transition between the 
ferromagnetic phase and the paramagnetic one. By in- 
serting it in (|43|) we obtain the free energy 

1 ~ \ if T >T C [ °> 

which shows a discontinuity in the second derivative at 
T = T Cl coherently with the continuity of the phase tran- 
sition. 

The picture of the phase transition is identical to the 
one of the mean-field Ising model apart non-substantial 
differences in the shape of the entropy. 



V. CONCLUDING REMARKS 

In this article we have introduced a classical spin model 
showing a first order Z2-symmetry breaking phase tran- 
sition. This model can be considered a drastic simplifi- 
cation of the mean-field Ising model pj} which conserves 
the minimal conditions sufficient to induce a phase tran- 
sition. These conditions are a double-well potential as a 
function of the magnetization to with the gap between 
the absolute minima proportional to the number of de- 
grees of freedom N. The difference between the model 
here introduced and the mean-field Ising model is the 
type of phase transition: first order in the former, con- 
tinuous in the latter. This is due to the discontinuous 
and continuous behavior of v(m) in the former and in 
the latter model, respectively. 

The double-wellness of the potential opposes to the 
concavity of the entropy of configuration space in deter- 
mining the shape of the free energy /(to, T) — v(m) — 
Ts(m). Indeed, the competition between these two 
shapes takes place giving rise to a double- well /(to) or a 
convex /(to) depending on the value of the temperature 
T: if T < T c then the shape of the potential wins and so 
the symmetry is broken, conversely if T > T c then the 
shape of the entropy wins and so the symmetry remains 
unbroken [l5l [T6j . 

This analysis has been possible due to the fact that 
the potential can be written as a function of to, more 
generally the order parameter. In this way the tendency 
to the symmetry breaking induced by the double- well po- 
tential and the tendency to leave unbroken the symmetry 
by the entropy of configuration space have been able to 
be separated. 

We think that this mechanism has good probabilities to 
be the general one acting whenever a symmetry breaking 
phase transition occurs, but we have to take into account 
that the potential of the model studied in this paper are 
infinite-range. This is not the most general case occur- 
ring in general physical systems, especially in classical 
physics. Nevertheless, there may be the way to general- 
ize the mechanism to finite-range systems. This is under 
study. 

The paper ends with the revisiting the solution of the 
Berlin-Kac model [H, [l9| (known also as the spherical 
model) in the mean-field version. The study reveals the 
same picture already found out in the mean-field Ising 
model, apart non-substantial differences in the shape of 
the entropy. 

If the potential is not expressible as a function of 
to, but always infinite-range there are evidences that a 
double-well potential with a gap proportional to N can 
be again at the base of a Z2-symmetry breaking phase 
transition, as suggested by the studies contained in [l|, Q . 
In that papers some models with first order and continu- 
ous phase transition are built by modeling a double-well 
potential generalized to iV-dim with gap proportional to 
N. 
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